MOMENTS OF THE TRANSMISSION EIGENVALUES, PROPER 
DELAY TIMES AND RANDOM MATRIX THEORY I 
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Abstract. We develop a method to compute the moments of the eigenvalue densities 
of matrices in the Gaussian, Laguerre and Jacobi ensembles for all the symmetry classes 
(3 £ {1, 2,4} and finite matrix dimension n. The moments of the Jacobi ensembles have 
a physical interpretation as the moments of the transmission eigenvalues of an electron 
through a quantum dot with chaotic dynamics. For the Laguerre ensemble we also evaluate 
the finite n negative moments. Physically, they correspond to the moments of the proper 
delay times, which are the eigenvalues of the Wigner-Smith matrix. Our formulae are well 
suited to an asymptotic analysis as n — > cx). 
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1. Introduction 



1.1. Background. Over the past twenty years Random Matrix Theory (RMT) has pro- 
vided a powerful tool to investigate quantum properties of electronic transport through 



ballistic cavities (quantum dots) [6|-|8| 



15, 31 



The purpose of this paper is to compute averages of the form 



(1) Mf' 



{k,n) = ^ I ■ ■ ■ I i^x^AYl^tiixj) Yl \xk - Xjfdxi- ■ -dxn 

I I \j=l J j=l l<j<k<n 



for finite n and k and for any value of /3 G {1,2,4}. Here S labels one of the Gaussian 
(G), Laguerre (L;,) or Jacobi (Ja,b) ensembles and the value of (3 corresponds to ensembles 
of real symmetric (/3 = 1), complex hermitian (/3 = 2) or quaternion self-dual matrices 
(/9 = 4). The function Wi3{x) is the weight of the ensemble: 

{^-13x^/2^ / = (—00,00), Gaussian ensembles, 

2,/3/2(fe+i)-ig-/3x/2^ J = [0,00), Laguerre ensembles, 

_ 3,)/3/2{a+i)-i^ J ^ j^pQ^i ensembles, 

where C is a normalization constant which may vary at each occurrence. The averages ([T]) 
for the Jacobi ensembles correspond to the moments of the transmission eigenvalues of the 
electric current through a ballistic cavity; the negative moments of the Laguerre ensembles 
are the moments of the density of the eigenvalues of the Wigner-Smith time-delay matrix. 

The physical dimensions of mesoscopic systems are such that the quantum nature of 
the electron becomes important and a classical treatment of its dynamics is not accurate 
anymore. Furthermore, at low temperature and voltage electron-electron interactions can 
be neglected; therefore, the electron scatters elastically inside the cavity, which is attached 
to two ideal leads connecting two reservoirs in equilibrium at zero temperature. If the 
leads support m and n quantum channels respectively, all the information on the electric 
transport is contained in the scattering matrix 



(3) S 



f' 

mxn 



tnxm 



The sub-blocks rmxm and tnxm are the reflection and transmission matrices through the 
left lead, while t^xn ''^'nxn those through the right lead. Without loss of generality, 
throughout this paper we shall assume that m > n. Since the scattering is elastic S is 
unitary. This is known as the Landauer-Biittiker formalism. 

The scattering matrix depends in a complicated way on macroscopic parameters like 
the energy of the electron and the shape of the cavity. If the classical dynamics inside 
the cavity is chaotic, then the fundamental assumption is that the electric current displays 
universal features; thus, it is natural to model the scattering matrix S with a random matrix 
drawn from one of Dyson's circular ensembles: the Circular Unitary Ensemble (CUE) when 
/3 = 2; the Circular Orthogonal Ensemble (COE) when /3 = 1; and the Circular Symplectic 
Ensemble (CSE) when /3 = 4. Let K denote a time reversal operator. If the dynamics is 
not time-reversal invariant then /3 = 2; if it is time-reversal invariant then /3 = 1 if = 1 
and /3 = 4 if K"^ = -1. 

In this paper we give a unified approach to compute the family of integrals ([T]) for all the 
(3 G {1,2,4} and give particular emphasis to those connected to statistics of the electric 
current. Our formulae are exact for finite matrix dimension. Since experiments can now 
be performed in quantum dots with a number of channels arbitrarily small |46| , recently 
there has been an increasing interest in computing finite n formulae 32|, |4J, |47|, l5J| . Some 
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of these i ntegrals have never been computed before, others are already available in the 
literature 
original. 



27, 30, 35, 44, 54 



In particular, most of the formulae for (3 = 1 and /3 = 4 are 
In §2] we will discuss in detail our results and specify which of the averages ([T]) 
are already known. 

Our formulae have distinctive advantages. Firstly, we can compute 'negative' moments 
for the Laguerre ensemble. Since the joint probability density function (j.p.d.f.) of the 



inverse delay times coincides with that of the Laguerre ensemble [16|, we obtain the mo- 
ments of their density. Secondly, for positive moments the sums in our formulae extend to 
the order of the moments k and not to the dimension of the matrices n. The sums that 
express the negative moments in the Laguerre ensemble run to n, but their limit as n ^ oo 
can be computed with little effort. As a consequence, although still relatively involved, our 
expressions are simpler and more manageable than those in the literature. Furthermore, 
our formulae provide a bridge between finite n results and their asymptotics. Indeed, in 
the second part of this work 41|, we compute the first three terms of the expansions as 



— 7- oo of the moments of the transmission eigenvalues and of the delay times. They agree 



with those recently obtained semiclassically |9Mll 



1.2. The Transmission Eigenvalues. The eigenvalues Ti, . . . of the matrix tt"^ are 
the transmission eigenvalues. The unitarity of S implies that the Ti, . . . ,T„ all lie in the 
interval [0, 1]. The dimensionless conductance at zero temperature is defined by 

(4) G:=Trtt+ = TrtY^ = Ti + --- + T„. 

Furthermore, if S belongs to one of Dyson's circular ensembles, then the j.p.d.f. of 
Ti, . . . ,Tn is 

1 " 

(5) P^^HTu...,T^) = -^ll^r n \Tk-T,f. 



C 

j=l ^<j<k<n 



The parameter a = |(m — n + — 1 measures the asymmetry of the quantum channels 



in the leads. Formula (EJ was first computed when m = n by Baranger and Mello [6| and 
by Jalabert et al f3l|] ; when m ^ n it was reported in this form by Beenakker [sl , where it 
was attributed to unpublished work by Brouwer (1994); i ts g eneral derivation appeared in 
the literature for the first time in an article by Forrester |23| . 

In a classic article Dyson [21j] classified complex many-body systems according to their 
fundamental symmetries and proved that they correspond to the random matrix ensembles 
labelled /3 G {1,2,4}. Zirnbauer [Hg} extended Dyson's classification scheme to Cartan's 
symmetric spaces and introduced new symmetry classes in Random Matrix Theory. Zirn- 
bauer also argued that these non-standard ensembles appear in the stochastic modelling 
of ballistic cavities in contact with a superconductor. Such mesoscopic systems are called 
Andreev quantum dots. In his PhD thesis Duefiez |2p| further generalized Zirnbauer's clas- 
sification. Furthermore, Altland and Zirnbauer j4| divided the symmetries of Andreev 
quantum dots into four fundamental classes. These ensembles are labelled by two integers 
{(3,6): as for Dyson's ensembles, (3 takes values in {1,2,4}; instead 6 G {—1,1,2}. The 
four classes are (1, —1), (2, —1), (4, 2) and (2, 1); they correspond to different combinations 
of time-reversal and spin-rotation symmetries. 

Our formalism applies to Andreev quantum dots too. Indeed, in a recent paper Dahlhaus 
et al 1^ computed the j.p.d.f. of the transmission eigenvalues. It is obtained by deforming 
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the right-hand side in equation (jS]): 

1 

(6) p^^''\T,,...,T^) = -l[T;'{l-T,f' n 



C 

j=l ^<j<k<n 



\n-T,\ 



As for the j.p.d.f. in a = ^ {m — n + 1) — 1. 

Equations ([5]) and are particular cases of the j.p.d.f. of the eigenvalues of matrices 
in the Jacobi ensembles, namely 



(7) rf!;,(xi,...,x„)=inxf(^-^^)-(i-x,o''/^(-^^^^ n 

j=l ^<j<k<n 

for < Xj < 1, (j = 1, . . . , n). We recover (j6]) by setting 

1 H — 1 — 1 and h = m — n. 



/3 V 2, 

The moments of the density of the transmission eigenvalues are defined by 
(9) te)-(Tr[(ttt)'=]) = Mj«(^, 



where a and b are given in equation ([8]). From a physical point of view, they are important 
because they are connected to the cumulants {(nj)) of the charge transmitted over a finite 
interval of time by the generating function |36j 

CO j oo / ^ \ ^ 

(10) E7(('^.)) = -E^(^^)(^^-i)'- 

j=l k=l 

(See also [l2|, appendix A.) For simplicity in this formula we have omitted the dependence 
on (/3, 6) and on the numbers of quantum channels m and n. The charge cumulants can be 
directly accessed in experiments [13] ■ Using the results in ^and the generating function 
in ffTOj) we can compute the cumulants to any given order. For example, the variance and 
skewness are given by 



:ila) {{k,)) 



nm(^ - 1 + n)(^ - 1 + m) 
- 1 + n + m){^ - 2 + n + m){^ - 1 + n + m) 



fllM IM = (n-m-^ + l)(n-m + i±^-l) 

^ ' {{k2)) (n + m + ^-3)(n + m + ^-l)" 

The special case 5 = of equations (ITTi) were computed by Savin et al js^] (see also fisl]). 
1.3. The Wigner-Smith Matrix. The Wigner-Smith time-delay matrix is defined as 

(12) Q = -^hS-'^. 

The individual eigenvalues ri, . . . , r„ of Q are called proper delay times, and their average 

(13) rw = -TrQ 

n 

is referred to as Wigner delay time. Here n is the total number of quantum channels in 
the leads. The Wigner delay time measures the extra time an electron spends in the cavity 
as a result of being scattered. If 5* belongs to one of the circular ensembles, then it was 
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shown by Brouwer et al IgJ that the j.p.d.f. of the inverses 'yj = Tj ^ {j = 1, . . . , n) of the 
proper delay times is 



n 

(14) P,{lu....ln) = ^\{lt"e-^^'^^'^ n 



C 

j=l l<j<k<n 



l7fc - Ij 



1/3 



where th is the Heisenberg time. In our context th = n. In a sequence of papers Savin 



and collaborators [48|, |49|, |5l| computed the probability distribution function of the proper 
delay times. 

The moments of the density of the proper delay times are defined by 
(15) (vl^'l) = -{TtQ') =n'-'Ml^J{-k,n), < n/3/2 + 1, 



n 

where in this case b = n— 1 + 2/13, and the right-hand side of (H^jl denotes negative integer 
moments of the Laguerre ensemble. 

The outline of the paper is the following: in §2 we present our main results; §3 is 
devoted to ensembles with unitary symmetry; in §4 we discuss the approach underlying the 
computations of the moments for ensembles with symplectic and orthogonal symmetries; 
§5 and §6 contain the proofs of the formulae for ensembles with symplectic and orthogonal 
symmetries, respectively. 

In the final stages of the preparation of this article, and after the results in this paper 
had been presented at two workshops j3 we received a preprint by Livan and Vivo js^] , in 
which some of our formulae were derived with a different method and approach. Their 
expressions are different, but equivalent to ours. 

Acknowledgements 
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2. Statement of Results 

2.1. The Moments of Transmission Eigenvalues and of the Proper Delay Times. 

Since the transmission eigenvalues are distributed with the j.p.d.f. of the Jacobi ensemble, 
their moments are the moments of the eigenvalue density of this ensemble for finite matrix 
dimension. A A-th order recurrence relation for the exact moments at /5 = 2 were first 



reported Ledoux [3J|. Explicit formulae were then obtained by Novaes [44| and by Vivo 
and Vivo js^l- Bai et al [Hj computed the leading order term of the asymptotic expansion 
for /3 = 1. 

From equation ( 11^ . computing the moments of the proper delay times is tantamount 
to calculating the negative moments of the eigenvalue density of the Laguerre ensemble 
for finite n. These negative moments have never been determined before, though positive 



integer moments were calculated by Hanlon et al [29[ and Haagerup and Thorbj0rnsen |28l |. 

It is worth reminding the reader that the moments of the proper delay times exist only 
for 

k<- + l, 

because for larger k the integral Ml^(— A;,n), with 6 = n— 1 + 2//3, diverges. 



^Random Matrix Theory and Its Applications I, MSRI, 13-17 September 2010, Berkeley, USA; VI 
Brunei Workshop on Random Matrix Theory, 17-18 December 2010, Brunei University, UK. 
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The general formulae for moments of the Jacobi and Laguerre ensembles are reported in 
the main sections of the paper. Throughout, the notation = r{n + k)/r{n) refers to 

the Pochhammer symbol; the binomial coefficient can take arbitrary complex arguments, 
i.e. 

, . ,'k\ T(k + 1) 

(16) ' * y ^ ) 



J J r{k-j + i)r{j + iy 

and is defined for negative integers by the limiting form 



(17) 



2.1.1. Broken Time Reversal ((5 = 2). 

Theorem 2.1. The moments of the transmission eigenvalues and of the proper delay times 
for f3 = 2 are 

fc-l min(j,n) / .x / . x 



and 



where 



(^/2 + m + n~2i+j + l){5/2 + m) (m) 



{5/2 + m + n-i) ^^^^^ [5 /2 + m + n - i + l) ^^^{5 /2 + n + l) {n + l) ^^^^ ' 



Remark 2.2. The moments for the Laguerre Unitary Ensemble can be defined even when 
k is complex and have the following particularly simple expression: 



(20) 

j=0 



(6 + n)(fc_j+i) 
+ n)(_,-) 



of which (fT9l) is a special case. If /c is a positive integer, the sum in (120|1 consists of at most 
k terms. 

Remark 2.3. The coefficients 

I fk\ f k 



(21) N{k,j) 

in formulae (|T8i) and (120|) appear frequently in enumerative combinatorics, where they are 
called Narayana numbers. 

2.1.2. Conserved Time Reversal with = — 1 ((5 = 4j. 
Theorem 2.4. The moments of the transmission eigenvalues are 

min(H,Lfe/2j)min{fc-2j,2n-2j) , . , . 

(22) (r,'«) = i(7^:S)- X: E (g(,4.)5i.(*:.,n.„) 
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The coefficient Sf^j{k,m,n) is 



{6/2 + 2n-z- 2j) (2m) {6/2 + 2m-l) 



(5/4 + n + 1/2) (m) (5/4 + m - 1/2) ^^.^.^ {2n - 2j + l) {n + l) ^^^.^ 

{6/2 + 2m + 2n- 4j) (5/2 + 2m + 2n - 2i - 4j + A;) 
(5/2 + 2m + 2n-t- 2j) ^^^^^ {6/2 + 2m + 2n-t- 4j) ^^^^^ " 

Furthermore, the moments of the proper delay times for /3 = 4 are given by 

(23) {v^^l)=n'-'Mi!l^^i-k,n), 
where the moments of the Laguerre Symplectic Ensemble are 

Mt^{k,n) = 2-'-'Mil{k,2n) 

(24) ^^:^^fk\f k ^^(26 + 2n)(fc_,_2,+l)(2n-^-2J + l)(,) 

+ 2j. 



j=l i=0 



The symbol [-J denotes the integer part. 

Remark 2.5. The order of the moments k in equations (l22l) and (l23l) is a positive integer. 
However, (12^ holds even when k is complex. As for /3 = 2, if A; is positive the sum in 
equation (124|) contains only A; terms. 

Remark 2.6. Although n is typically an integer, as it denotes the dimension of a matrix, 
the expressions on the right-hand sides of equations fl22l) and fl2^ are well defined for any 
half-integer n. It is useful to generalize it, because the evaluation of the moments for /3 = 1 
requires moments for /3 = 4 computed at half-integer n. 

2.1.3. Conserved Time Reversal with = 1 ((5 = 1). For simplicity we assume that the 
outgoing lead supports an even number of open channels. 

Theorem 2.7. The moments of the transmission eigenvalues are 

min(?i/2— l,fc) , ^ 

(25) (ri^;S,) = 2 + E (((:) I^ik, m. „) + 4,1 

where 

(26a) /f'^)(fc,m,.)=4^ (5 + m + n - 4, + 2A:)(i(5 + m + l))(,_,)(lm)(,_,) 



and 



(26b) <j^i^ = J2 



{6 + m + n- 2j)(2fc+i)(|(5 + n + 2))(_,)(i(l + n))(_,) 
2^+2r(|(5 + m - n + 2j + l))r(i(5 + m + 2)) 



r(5 + 1 + m + J + k)T{j + k + l- n)r(i(m - n + 1 + 2j)) 
r(i(5 + n + 2))r(m -n + k + j)T{j + k) 



X 



r(f)r(| + i)r(f; 



r/ie moments of the proper delay times are 

(1) \ _ ^fc-lol-fc /\^(4) 



^^i::; = ^'^-^2-'=<;U/,(-fc, (^-i)/2) 

^^^^ ^ '^-^ "y^' ^2A: + 2j - 1\ (n + l/2)(_,_,) , 

^2; 2j J 2(1(1 + .))(_/" 
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where 

,L ^ r(n) ^ r(fc + j + n)r(l + n-fc-j) , 

^ ^ r(n/2)r(2n)^ r(n/2 + i-j)r(fc + j + i) " 

Remark 2.8. Due to the term T{j + k + 1 — n) m the denominator of fl26bp . the (pi^s 
are identically zero for any n > 2k. By Stirling's formula, the 0^1^ ^'s in fl28|) decay expo- 
nentially fast as n — > oo. Therefore, neither of these terms contribute to the asymptotics 
of the moments as n — )■ oo at any finite algebraic order. 

It is straightforward to compute the limit as n — )■ oo of the formulae in this section. 
They differ fundamentally from most of the known exact results in the literature, whose 
asymptotic limit often involves many cancellations, which means that even the leading order 
term may be out of reach. This difficulty is discussed in some detail by Krattenthaler 



where a solution is presented for /3 = 2 (see also [13]). 

Indeed, it is a simple exercise using our exact results to show that 

(29) lim i(7;,f,S)=(l + ag(*-l)Q(-l)V«. 

j=o ^ / 

where f is the variable f = , ""^no , which remains finite as n,m —)■ oo, and C, = ^tt(^'^) 
is the j-t/i Catalan number. This formula agrees with the semiclassical computation of 
Berkolaiko et al [9i] Furthermore, for the proper delay times we have 



\2K 



(30) i-?iK)^lgG)G-:)^ 

which is the k-th Schroder number (note the appearance of the Narayana numbers (12T|) ). 
This limit was computed semi classically by Berkolaiko and Kuipers (lO|, and can also be 
obtained from the Marcenko-Pastur distribution [s^ (see, e.g., [lo|, [l6i]). It is a simple 
consequence of (IT^ too. 

Equation fl29|) was first computed using RMT by Novaes jisj (se e also jl]), while f l30|) and 
fl29l) were recently rederived through combinatorial techniques j45| . Our exact results allow 
a simple derivation of these facts, while also consenting the investigation of /3-dependent 
subleading corrections. We address these issues more thoroughly in the second part of this 
work j4l|, where we show that the first two subleading terms in the asymptotic expan- 
sions of the previous theorems agree with those obtained semiclassically by Berkolaiko and 
Kuipers jllj . 

2.2. The Gaussian Ensembles. Our techniques apply equally well to the Gaussian en- 
sembles. Recursion formulae for the finite n moments of the density of the eigenvalues 
were derived by Harer and Zagier 130] for the Gaussian Unitary Ensemble (GUE), while 
Goulden and Jackson [22] derived explicit formulae for both the Gaussian Orthogonal En- 
semble (GOE) and the GUE, while the GUE moment generating function was computed 
by Haagerup and Thorbj0rnsen [2^. More recently, recursion formulae were obtained by 
Ledoux l35| for the GOE and GSE. 



Theorem 2.9. The moments of the eigenvalue density for the GUE are 
(2). , . 2"r(n/2 + l)r(n/2)™"^;^~'''V'^\ A + . 

^3^) ^s\2k,n)^ v^(k+i)r(4 g G)G+i)(^/'-^'^^-^/^^ 
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for even n, and 

/or oc?(i For the GSE we /iave 
(33) 

M^\2k,n) = 2-''-^M^(^\2k,2n) 

Let n be even. Then, the moments for the GOE are 
M^^^ {2k, n) = M^^^ {2k, n - 1) 



min(f-l,fc)min(fc,f-l-j) ..x , , x /„ . .n 

'k\f k \[2-^-J)(k+i/2) 



E E 

j=l i=0 



Ij\i+Jj (f-j)(l/2) 



For n <2k the quantity is given by 



fc-n/2n/2-l /n-l\22^2!±:i)i 



(2A:)! "^^^ (n/2 -.-!)! 
r(n/2) ^ ^ (^■_,)!(A:-j-)!4^- 



If n > 2k we have 



(34) 



^.^^ 2'=-3.(2j)!(A;-j)! 



3. Unitary Ensembles 

We shall now compute the moments of the eigenvalues densities for the Jacobi, Laguerre 
and Gaussian ensembles when /3 = 2. For brevity we shall refer to these ensembles with 
the usual notation JUE, LUE and GUE. Except for the GUE, our expressions are valid for 
complex k. Theorem 12. II and equations (IHTl) and ([32]) of theorem 12.91 are corollaries of the 
results of this section. 

For all the ensembles and symmetry classes that we consider the j.p. d.f. of the eigenvalues 
has the form 

1 " 

(35) pf\xi,. . . ,Xn) = -^Wwp{Xj) JJ \Xk-Xjf. 

The marginal probabilities are obtained by subsequent integrations of the right-hand side 
of (1351) : furthermore, since it is invariant under permutations of its arguments, it is ir- 
relevant which variables are integrated over. Therefore, the probability density of the 
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eigenvalues is obtained by integrating out all but one variable. It follows that 
(TrX'^) = j--- j{x\ + --- + xl)pf\xi,...,Xn)dxi---dxn 

x''pi3{x)dx, 




(36) 



where P/s^x) is the eigenvalue density normalized to n and / is the support of W/3(a;). 

We develop effective techniques to compute the integral in the right-hand side of equa- 
tion f l5B]) using ideas first introduced by Haagerup and Thorbj0rnsen [28|] for /3 = 2 and by 
Adler et al for /3 = 1,4. 

When /3 = 2 the density of the eigenvalues takes a particularly simple form (see, e.g., (25| . 
§5.1) 

(37) p2{x) = 

where the Pj(a;)'s are orthogonal polynomials associated with the weight W2{x) and j = 
0, 1, . . . denotes their degree. In other words, we have 

(38) j W2{x)Pj{x)Pk{x)dx = hj6jk, j,k = 0,1,... . 

The system of orthogonal polynomials {Pj{x)}'jLQ is unique up to multiplicative constants 
kj, which we can take to be the coefficient of the monomial of highest degree. Orthogonal 
polynomials satisfy a recurrence relation of the form 

(39) P,+,{x) = {a,+x(3,)P,{x)-^,P,,^ix), j = 0, 1, . . . 

where for convention P_i(x) = 0. For the classical orthogonal polynomials the constants 
hj, kj, aj, Pj and 7j are tabulated in many books on special functions (see, e.g., [l|). A 
consequence of fl39|) is 



(40) p2ix) = W2ix) (P^(a;)P„_i(a;) - P„(x)i^:_i(x)) , 

which is a limiting case of the Christoffel-Darboux formula. (For the proofs of formulae (15^ 
and (gni) see, e.g., §3.2). 

In the rest of this article we shall assume that kj = 1, for j = 0,1, . . . . In other words, we 
only consider monic orthogonal polynomials. In order to distinguish them from the way the 
classical polynomials are conventionally defined in the literature, we shall use the notation 
Tinix), C^{x) and V^'^^i^x) for the Hermite, Laguerre and Jacobi polynomials respectively. 
We shall denote the generic monic polynomial by p{x). We tabulate the orthogonality 
constants hj for the monic classical polynomials in appendix A. We shall also need the 
following differential equations (see [1], §22.6) 



(41) 



'U'l^x) - 2xU'j{x) + 2jHj{x) = 0, 
xCj{x)" + (6 + 1 - x)C]{x)' + jC){x) = 0, 
x(l - x)V;'\xy + (6 + 1 - (a + 6 + 2)x)V;'\x) 
[+j{a + b + j + l)V;^\x)=0. 



Haagerup and Thorbj0rnsen [28| computed the moment generating function 
(42) M{t) = [ p2(a;)e-*^rfx 
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in terms of hj^ergeometric functions for the GUE and LUE. They combined the differential 
equations fllT]) with fHU]) to obtain 

d j -D^e~''^'Hn{x)'Hn-i{x), Hermite, 

ax 1— iJ^x^e ''£„(x)£„_^(a;), Laguerre, 

where 

(44) Z5H^___ and = (r(6 + n)r(n))-i. 

Furthermore, /(x) = 1 for the Hermite polynomials, while f{x) = x for the Laguerre ones. 
We shall use similar ideas to compute the moments (Q for /3 = 2. 

First we need the analogue of the identities 0431) for the Jacobi polynomials. 

Lemma 3.1. Let P2{x) be the mean eigenvalue density for the JUE. We have the following 
differential identity: 

(45) ^ (x(l - x)p2ix)) = -Dix\l - xrv:i'\x)V:'\ix), 
where 

J _ T{a + b + 2n + l)T{a + b + 2n - 1) 
^ ^ " ~ r(a + n)T{b + n)T{a + b + n)T{n) " 

Proof. The normalization coefficient hn-i associated with the polynomials Vn-i{x) is 

T{a + n)T{b + n)T{n)T{a + b + n) 



(47) h 



n-l 



T{a + b + 2n)T{a + b + 2n-l) 



Inserting ( H7|) into the representation (HOl) and using the differential equation in (HTj) we 
obtain 

(48) x(l - x) (P^) ' + (1 + 6 - (a + 6 + 2)x) (P^) = -Div:\x)V::\{x), 

\W2{X} J \W2{x)J 

where is given in fHB]) . Finally, since the weight associated with the Jacobi polynomials 

is ^2(3;) = x^(l — x)" we arrive at 

= x\l-xy (i{b+l)il-x) -x{a + l))^^ + xil-x' '^^2(a;)^ 



W2{x) \W2{x) 

-Dix\i-xrv:'\x)v:'\ix). 



□ 



Remark 3.2. For our purposes it is not helpful to compute the moment generating func- 
tion Although in principle one can employ a type of fractional calculus to extract 
more general types of moments from fH2|) . we will see in the following that moments for 
general k are directly accessible with our method. 
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3.1. Jacobi Unitary Ensemble. Lemma [3.11 allows us to compute the difference of the 
moments. Then, the moments themselves can be computed by adding all the differences. 
Finally, equation f[T5]) of theorem 12. II is obtained by setting a = 6/2 and b = m — n. 
Let us define 

(49) AMf = M^l\ik,n) - Mf\{k + l,n). 



n,a,b 



Proposition 3.3. We have 

(50) ^<(*.") = ^EG)(,!i)f' 

where 

(51) f/"'"''' = (Q + ^ + 2^ - 2j + + + & + n)(fc_j+i)(a + n-j + l)y){b + n)(^k-j+i) 

(a + 5 + 2n- j)(fc+2)(a + 6 + 2n-j + l)(fc)(n + l)(_j) 

If k is a positive integer, equation (l50l) reduces to 

(52) AMj:i(..„) = i g OC ^Of't- 

Proof. Integrating by parts using equation leads to 

(53) x^(l - a;)p2(a;)c/x = ^ f a;'=+'(l - xfV^J{x)V^^\{x)dx 
Jo k, Jq 

Consider the identity 

n 

(54) P^'^(x) = 5^C^'='"P;''+'=(x), 

i=o 

where 

,fc,n A\ (a + n + 1 -j)Q)(a + 6 + 2n + l)(_j) 



f55) C 

Vj7 (« + ^ + 2ri-2j + 2 + A;)(j)(n + l)(_j) 

are the connection coefficients. Inserting this formula into (!53|l and evaluating the integrals 
using orthogonality leads to 

(56) AMjJ) (fc,n) = C^^r'/^^'Y- 

j=0 

Substituting the appropriate coefficients (see appendix A) gives immediately (l50i) . 

When is a positive integer, the terms with j > k vanish because (^) = if j > /c, 
leading immedietely to f l52|) . □ 

Corollary 3.4. T/ie integer moments of the level density for the JUE are 

(57) M«(*.„) = Mj:)(l,„)-X:i E 

where the first moment is 

(58) Mf'd.n- + 



Tn,a,b 
hi ■ 



a + 6 + 2n 

The first moment M] (1,"^,) is an Aomoto integral. For its evaluation see, e.g.^ [39 
il7.3. 
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3.2. Laguerre Unitary Ensemble. Since tlie moments of tlie Wigner-Smitli matrix f[T^ 
require the computation of the integral fl5^ for A; < 0, we shall present formulae for the 
moments of the LUE for general complex k. 

Proposition 3.5. Suppose that neither b + n nor b + k are negative integers. Then one 
has 



k \{b + n)(^k-j+i) 



(59) 

j=0 

Proof. Integrating by parts the second equation in fj^3l) gives 



(60) / x^p2{x)dx = -^ / x''+V"£^(x)£^_i(a;)rfa;. 



For the Laguerre polynomials the connection formula is [52 



(61) 4(a:) = X:C^'"4^^x), where C^'" = (^) (r. + 1)(-,). 
Inserting formula fl6T|) into fl60|) gives 

(62) Mif (A;,n) = ^ J]cf C^r^Cj, 

j=0 

where we evaluated the integrals using orthogonality. Using the appropriate connection 
coefficients and normalisation constants completes the proof. □ 

If A; is a positive integer, the binomial coefficient = if j > fc, leaving only a sum 
with k terms. Negative moments are obtained simply by using the identity 



Corollary 3.6. Let k be a positive integer, then 

min(?i,fc) 



(63) M<:'(*,„) = i g (*)(^,*^ 

Furthermore, if k < n + 1 we have 



k \ (6 + n)(fc_j+i) 



(n + r 



i-i) 



Equation f|T9l) is a particular case of formula flMl) . where b = n and the scaling introduced 
by the Heisenberg time tu = n has been taken into account. 

Remark 3.7. The appearance of the Narayana coefficients in (l63l) anticipates the fact that 
its leading order term as n — )■ oo is the k-th moment of the Marcenko-Pastur law 38 . 
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3.3. Gaussian Unitary Ensemble. In this section we give the proof of equation ( 13 ip of 
theorem 12.91 The approach is the same as for the JUE and LUE. The proof for n odd is 
very similar and we omit the details. 

Integrating by parts the first formula in (H3|) gives 

/oo poo 
x^''p2{x)dx = / x^^^^e"'^nn{x)'Hn-i{x)dx. 

oo 2k -\- I J -oo 

The integral f lB^ can be evaluated using the Laguerre polynomials since 

(66) 'Hn(x) = £^/2^(x^) and Un-iix) = x£^^2-i(^^)- 
A change of variables then leads to 

jjU poo 

(67) Mg\2k,n) = x'+'/'e-^Cl!;,_,{x)C:ji\x)dx. 

This integral is of the same type as that one in the right-hand side of equation (l60i) and 
can be computed in the same way. 



4. Symplectic and Orthogonal Symmetries 

Very few non-perturbative results are available for the moments of the densities of the 
eigenvalues for /3 = 1 and (3 = 4. Our goal here is to develop a novel approach that allows 
us to compute these moments for all the ensembles associated with the weights ([2]). 

There are two possible ways of tackling this problem: the first is through the Selberg 
integral; the other one is a direct computation of the moments (l36i) . The Selberg integral 
was very effective in computing the moments of the transmission eigenvalues for /3 = 2 



when /3 = 1 it does not seem to produce explicit formulae 3^. It cannot be applied to 
/5 = 4. 

Following an approach of Dyson j,22|], Mehta and Mahoux [4C!j expressed the densities 
for (3 = 1 and /3 = 4 in terms of skew-orthogonal polynomials. Since then several articles 



have attempted to improve their formulae [2|, 126|, 142, 153|, 155(. Tracy and Widom |53|^ 
and Widom [ss!] succeeded to write such densities as sums of P2{x) plus correction terms 
involving orthogonal polynomials. Building on the work of Adler and van Moerbeke j^, 
Adler et al [2] obtained integral representations of the correction terms. 

Equation fl36l) presents one major challenge: for finite n it is a complicated sum involving 
all the orthogonal polynomials up to n — 1. Further integration would lead to cumbersome 
formulae whose asymptotics cannot be easily extracted. Our method relies on using the 
coefficients flM|) and fl6Tl) to expand the orthogonal polynomials in a convenient basis, 
within which they are orthogonal with respect to the perturbed weight x^W2{x). As when 
P = 2, this allows us to obtain positive moments involving sums that run to the order of 
the moments and not to the dimension of the ensemble. Another interesting feature of our 
results is that we are able to express the moments at /3 = 1 in terms of the moments at 
/3 = 4 plus a fairly simple correction term. Like for unitary ensembles, our formulae are 
sum of ratios of Gamma functions which may be studied in the limit n — )■ oo. 

Since our approach is based on the results in by Adler et al [2], we will discuss their 
formalism in detail. For /3 = 1 and /3 = 4 a special role is played by the skew-orthogonal 
polynomials. Recall that an inner product {A, B) is referred to as skew if (A, B) = —{B, A). 
A sequence of monic polynomials {gj(x)}°?,o called skew-orthogonal with respect to 
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(A, 5) if 

(68a) {Q2m,<l2n+l) = — {Q2n+1 , Q2rn) = ^m^m,n 

(68b) {q2m,q2n) = {q2m+l , q2n+l) = 0. 

Let us introduce the potential V{x) by defining 

(69) W2{x) = e-^^^'\ 

where W2{x) is the weight function of the associated unitary ensemble. We also assume 
that 

(70) 2V'ix) ^^^^ 



(72a) {A, B)i = l I e-2^''(^) {A{x)B'{x) - B{x)A\x)) dx 



fix) 

is a rational function of x and take f{x) to be a monic polynomial. Now define the modified 
potentials: 

(71a) V,{x) = V{x) + hogf{x), (3 = 1, 

(71b) V,{x) = Vix)-hogfix), 13 = A. 

Let us introduce the inner products 

and 

(72b) {A,B)i = ]^ j J e-^^^''^-^'^yhgn{y - x)A{x)B{y)dxdy. 

Associated to these inner products are two systems of monic skew-orthogonal polynomials 
{qf\x)}^Q and {qf\x)}'^Q. We shall denote their skew-norms as defined in ([2HD by fj"^^ 

and respectively. The tilde notation indicates that the weight has been perturbed by 
the transformations fl71al) and fl71bp . 

Because the skew-orthogonality relations fl68|l are invariant under the transformation 
q2m+i — ^ 92m+i + a2m9^2m for any a2n e C, a system of skew-orthogonal polynomials is 
not uniquely defined. However, they can be expressed in terms of the monic polynomials 
orthogonal with respect to ^2(0;): 



(73a) P2j+i{,x) = qt-^i{x), p2j{x) = qif{x) - ^^^qi^^ 



x 



(73b) qif{x) = P2j{x), qif+i{x) = P2j+i{x) - ^^p2j-i(x). 

For the classical orthogonal polynomials the constants in these equations are given by 

(74) Cn = K+lK^n, 

where 

{1, Hermite, 

i, Laguerre, 

i(2n + a + 6 + 2), Jacobi. 

We point out that the numbers appearing in the differential identities may be 
expressed in terms of 7^ via = 2^^-! for all the three ensembles. For each ensemble, 
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we denote the mean eigenvalue density by Pi3{x), where the tilde indicates the ensemble 
average (15^ defined by the weight e~'^'^^^^ for /3 = 1 or e~^^*(^) for /3 = 4. 

Remark 4.1. Before we proceed it is worth noting that the weights e"^^*^^'-* turn out to be 
exactly equal to the weights wi{x) in equation ([2]). The eigenvalue densities p4(x), however, 
correspond to the weights e~^^''(^\ which are 

{e~^^, Hermite, 
x^'+^e-^ Laguerre, 
Jacobi. 

These are not quite the same as the weights in ([2]) for /3 = 4. We must make the substitution 
(a, b) — )■ (2a, 26) for the Jacobi ensemble and b ^ 2b for the Laguerre ensemble. In 
addition there is a missing factor of 2 in the exponentials which we take into account 
by multiplying our final results for the moments by the appropriate power of 2. This 
discrepancy arises because the symplectic ensembles are sets of self-dual n x n quaternion 
matrices. Their representation in terms of complex matrices leads to Kramer's degeneracy, 
which is responsible for the normalisations in fl75]) . Without loss of generality we shall still 
use the notation ^4(0;). 

Let us introduce the e-transform of a suitable function f{t) by 
(76) e[f{t)]{x) = ^ l^sgnix- t)f{t)dt. 

In the following, / = (mi, 7712) will denote the interval of orthogonality. We have the 
following formulae 

/•m2 

(77a) p^ix) = -P2{x)n^2n " 7;l2n-ie-^'^''^ P2n{x) / 6"^^ Wp2n-1 (i)o?^ 



and (for n even) 

(77b) pi{x) = P2{x)n^n-i + 7n-2e-^^ Vn-1 (x)e [p„_2 (t)e-^i (x) . 

For convenience we shall alter the representations (I77al) and (177b|) into a form which is 
more suitable for the evaluation of the integrals 0361) . The following proposition allows us 
to expand the integrals in fl77al) and fl77bl) in terms of monic orthogonal polynomials. 

Proposition 4.2. Let {pj{x)}'^Q be the system of monic polynomials orthogonal with re- 
spect to W2{x). We have the following identities: 

n 

(78) e [e-^^Wp2n+2(t)] (x) = -e'^^^^) ^ e«p2,+i(x) + [e'^^W] (x) 

i=o 

and 



(79) / e-^^Wp2.+i(t)rft = -e-^*(^) J] e'flp2j{x), 



where 



\ 71, Tl 

^0) efi = ^r\'-^, e« = ^n^' r^« = n^^%^- 

C2n C2i ' C2,+l . C2j+l „C2j+lh2j 
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Proof. We begin from the differential identities [2| 
d 



(81a) f (e-^^'M%)) = T^e-^^(^Wi(x), 

dx \ J /l2n+l 

(81b) ^ (e-^^(-)gtVi(x)) = e-^^(^) (^P2n(x) - ^P2„-.2(x)) . 

We also need e'^^'y^^^ = e~^4(m2) _ which can be easily checked from fl75]) . 
We first derive ( 17^ . Integrating equation (18 lap between x and m2 gives 

(82) / e-^^«p2n+i(t)rft = -^e-^^^%!£{x) = -e'^^^^) J] e(.%(a;), 

where the last equality was obtained by iteratively solving equation f l73ap for qini^)- 
In order to derive (ITS]) , we start by integrating equation (181bp between x and 

(83) r e-^^(S2„+2(t)rft = '-^^ r e-^^^%,^m - ^e-^^^^^^'Ux). 

Jx C2„+l/i2n Jx C2n+1 

Integrating (ISlbp between mi and x and subtracting the result from (|83|) gives an equation 
for the e-transform 

(84) e [e-^^Wj,2.+2(t)] = "-^^e [e-^'^^'Wit)] - ^e-^^(^W.+i(x), 

C2n+l"'2n C2n+1 

where we used that qif+i — P2j+i{x). Iterating this equation n times leads to (ITHj) . □ 

Remark 4.3. The coefficients e^^l, e^j^l and rjn^ are tabulated in appendix A for each 
ensemble. 



Corollary 4.4. We have the following representations for the eigenvalue densities: 

(85a) p^ix) = -p2{x)n^2n " -72^-16"^^^"^^ ^ e''!^l_^P2jix)p2n{x) , 

j=0 
n/2-2 

(85b) pi(x) = P2(a;)n-fn-l - 7n-2e"^^*^^^ <^^ll/2~2P2j+li^)Pn~l{x) 

j=0 



+ 7n-2e-^^(^'Vn-i(a;)r/i/L2e (x). 

Proof. Substituting the integration identities (1791) and (1781) into (I77al) and (]77bp respec- 
tively, and using Vi{t) + V^it) = 2V{t), gives (l85al) and (l85bl) . □ 

We will see in ^ and ^ that formulae (185|) are particularly suited to our purposes. A 
key feature of these representations is that they are expressed solely in terms of the weight 
function e~'^'^^^^ and the corresponding monic orthogonal polynomials. For the orthogonal 
ensembles, there is an additional term involving the e-transform of the weight e~^^*^*\ which 
is related to the error function, incomplete gamma function or incomplete beta function 
depending on the ensemble in question. We shall compute the moments for (3 = 1 and 
/9 = 4 by combining the representations (1851) with a variant of the technique used for 
unitary ensembles. 
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5. Symplectic Ensembles 
The purpose of this section is to compute the integrals 

(86) M^^\k,n) = j^x''pi{x)dx 

for each ensemble S defined by the weights ([2]). 

Inserting the representation fl85al) into equation (l86l) leads to two integrals: the first one 
contains the mean eigenvalue density of a unitary ensemble, which was computed in ^ 
the second one involves orthogonal polynomials. More explicitly it is given by 

(87) Ss{k, n) = 72„_i- ^ e^.J.i / x''e-'''^%,,{x)p2nix)dx. 

j=0 -^^ 

We know how to evaluate these integrals as they are exactly of the type that appeared in 
equations (153|) . (160|) and (l65|) . We write the polynomials Pn{x) in a basis which is orthogonal 
with respect to the perturbed weight x^e"^^*^"^-*; then, we can use the orthogonality of the 
polynomials to write the integral in f E7|l as a single sum involving the connection coefficients 
(jSlj) and fl^T]) . Eventually, the moments for /3 = 4 become 

(88) MP{k,n) = Mf\k,2n) - S£{k,n). 

As in ^ the tilde notation indicates quantities that differ from the integrals ([T]) by a 
factor discussed in remark 14.11 



5.1. Jacobi Symplectic Ensemble. We now compute Sj^^{k,n) for complex k. Equa- 
tion fl2^ is then obtained by restricting /c to be a positive integer and setting a = 5/4 — 1/2 
and b = m — n. 



Proposition 5.1. We have 

n 2n-2j 



j=i 1=0 V 



where the coefficient S°''j{k,n) is given by 
(90) 

ga,b^f^ . ^ 2^^~\2a + 2n-i-2j + 1)(,)(2& + 2n)(fc_,_2,+i)(2a + 2b + 2n)(k-i-2j+i) 
' (2n - 2j + l)(_i)(n + l)(-i)(a + n + l)(-j)(& + ^)(i-j)(« + ^ + ^)(i-i) 

(2a + 26 + 4n - 4j + l)(2a + 26 + 4n - 2i - 4j + A; + 1) 



X 



{2a + 2b + An-i- 2j + l)(i+fe)(2a + 26 + 4n - i - 4j + l)(i+fc) ' 
Proof. By f l87|) we have 

(91) 5j„,,(A:, n) = ^ 5^ eJ5_, / x^+^1 - xrV-,;,\x)V^,;{x)dx. 
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Inserting the connection formula flM|) into the integrand leads to 

n-l 2j 2n „i 



j=0 1=0 p=0 

n-l 2j 2n 



(Qn\ _ \P „(4) 'Sp/^fc,2i V^^fc,2n^a,6+fc^ 

ly^J - 2^ Z^^i Z^^P "'2j-i (>2n-p,2j-i 

j=0 i=0 p=0 

n 2n—2j 

_ 72n-l \ ^ \ ^ (4) ia,b+k pk,2n-2j pk,2n 

~ 2 / ^ / J ^n-j,n-l'^2n-2j-i^i i+2j- 
j=l i=0 

To obtain the above expression we have applied the orthogonality of the Jacobi polynomials 
and then rearranged the indices in the sum. Substituting the coefficients h'^^^^ and 

e^^\ (see appendix A) and replacing (a, h) — )■ (2a, 26) completes the proof. □ 

Remark 5.2. The complexity of the expression (I^Uj) is mainly due to formula flM|) . which 
relates Jacobi polynomials of different weights. The Laguerre ensemble is slightly simpler, 
because the associated coefficients flM]) and normalizations hj are more concise. 

5.2. Laguerre Symplectic Ensemble. 

Proposition 5.3. Suppose that neither 2b + k nor b + n are negative integers. Then, we 
have 

N 2n-2i , ^ + 2n)(fc_i_2,+i)(2n - t - 2j + l)(i) 



Proof. From flHTj) we have 



(94) ^L,(A;, n) = ^ ^ ej^.i / x'^'e-£l{x)C',^{x)dx. 

Proceeding as in the proof of proposition 15.11 we obtain 

n-l 2j 2n „oo 

Suik, -) = ^ E eJt, Yl E C.''" / ^'^'e-^i^ll%{x)C\f_,{x)dx 

j=0 i=0 p=0 -^^ 

n-l 2j 2n 

j=0 i=0 p=0 

n 271— 2j 

_ 72n-l \^ (4) ^k,2n-2jpk,2n,b+k 

~ 2 / ^ / ^ ^n-j,n-l^i ^i+2j"'2n-2j-i 
j=l i=0 

By replacing b ^ 2b and multiplying both sides of this equation by 2~^ gives the statement 
of the proposition. □ 

Both the propositions 15.11 and 15.31 hold for complex values of k, except where the inte- 
grals ([T]) diverge. In particular, the moments of the proper delay times ([T^ are expressed 
in terms of negative moments of the Laguerre ensemble and can be obtained from 
using the identity 

(-k\ ( -k\ (k + 3-l\/k + , + 2] - r 
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and setting 6 = n + 1 in (jHSD- Thus, we arrive at the following: 
Corollary 5.4. Let k be a positive integer. We have 

min( [nj , [fc/2j ) min{2n~2j,k-2j) 



(98) 



j=l j=0 

(26 + 2n)(A._i_2i+i)(2n - i - 2 j + 



Furthermore, if k < 2n + 1 

5,,(-A:,n) = 5^5^ (, J (, A: - 1 J 

(99) i=i i=o ^ / ^ ^ 

(26 + 2r2)(,fc_i_2j+i)(2n - i - 2 j + 
X 2-'=-2.+2(^ + i)(_^.^(^ + 5)(^__^.^ • 

Remark 5.5. The combinatorial aspect of the sums in this corollary arises directly from 
the connection coefficients Cj and Cj+2j appearing in (jHSD and (1^ . leading to the binomial 
coefficients. Due to these binomial coefficients, the sums fl97|l only go up to k and their 
complexity does not increase as n grows. This therefore permits an investigation of the 
asymptotics. The — )■ oo analysis of fl99l) leads to certain infinite sums which also turn 
out to be tractable. A similar remark holds for the GSE. 

5.3. Gaussian Symplectic Ensemble. In the Gaussian case, only the even moments are 
different from zero. For simplicity we only state the results for integer k. 

Proposition 5.6. We have 

min(n,A;) min(n—j,k~j) 

I \ I \ / 

" ' ^ i(fc-l/2) 

j=l i=0 \ / \ ^J/ 

Proof. The integrals f l87|) give 

(101) ^G(2fc,n) = I^^efi„i / x^^e'^"n2n{x)n2,{x)dx. 
Applying formula (l66l) leads to 

rt— 1 „oo 

(102) ~So{2Kn) = ^ J^eJ-ti / x''e~'^' C-'^\x-)Cj'^\x^)dx. 
Changing variables and inserting the connection formula (16T!) results in 

n-l j n „oo 

j=0 i=0 p=0 -^^ 

Using the orthogonality of the Laguerre polynomials gives the double sum 

n n—j 

(103) Sa{2k, n) = ^-^J2Il ^.n-^'^'C^^i^r 

j=l i=0 



■P^ I 1\-P/ \ min(n,A:) min(n-j,fc~j) . , x 



20 



It is worth emphasising that the coefficients C^'^ ■' , C^^^ and hn are those of the Laguerre 
polynomials, while the coefficient e^^}_j^_i is related to the Hermite polynomials. Finally, 
multiplying equation fllOSp by 2~*^, as discussed in remark Kl\ completes the proof. □ 

6. Orthogonal Ensembles 

In this section we compute the moments for /3 = 1. Theorem 12.71 and equation (IMj) of 
Theorem 12.91 are corollaries of the results we prove here. For simplicity we assume that n 
is an even integer. 

The main task is to compute the integral 



(104) M^/'{k,n) = J x''pi{x)dx. 

When (3 = 1 the density pi{x) coincides with pi(x), so the integrals fll04p coincide with 
the averages ([T]). 

Substituting the representation (]85b[) into (11041) . we are left to compute three integrals: 
the ffist one gives the moments of the corresponding unitary ensemble; the second one is 
closely related to the quantity Ss{k,n) discussed in ^ namely 

n— 1 „ 

(105) 0£{k,n) = 7n-2X;4w2-2 / a;'^e-2^(^)P2,+i(x)P„-i(x)rfx; 

j=o -^^ 

the last one arises from the e-transform in equation fl85bl) . i.e. 

.(1) 



(106) h{k,n) = 7„_2<y2_2y^a;'=P„-i(a;)e-^^(^)e [e-^^^j {x)dx. 

Therefore, the moments for (3 = 1 may be expressed as 

(107) M^^\k,n) = Mf\k,n - 1) - Os{k,n) + Ie{k,n) 
In this section we focus on the integrals Os{k,n) and Is{k,n). 

6.1. A Duality between (3 = 1 and (3 = 4. We ffist discuss a remarkable duality between 
the quantities {k,n — 1) — 0£{k,n) and the moments of the symplectic ensembles 
M^'^\k,n), where n now can assume half integer values. Such moments are well defined 
(see equation (IHHD and remark 12. 6p . Similar dualities have appeared in the literature 
before (Til, M, SJ 



Lemma 6.1. Let n be an even integer. We have the following dualities: 

f^ik,n)=2'+''M^'^ 
r(i) tu ^\ _ o/i#(4) 



(108a) Mll\k,n) = 2i+X:;,(^, - l)/2) + Kik,n), 



(108b) K:,ik,n) = 2M];;^^^^(fc, {n - l)/2) + /j„,,(A;,n). 

Proof. Firstly, by equation f l88|) we observe that 



(109) 2'^'Mll^ik, {n - l)/2) = M^^^(A:,n - 1) - 25l,(A:, {n - l)/2). 



>i+'=<;^(fc,(n-l)/2) = Mif( 
Thus, it is sufficient to check that 2S'l^ (/c, (n — l)/2) = Oi^^{k,n). 
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A direct computation shows that 

j=o Jo 

n/2-2 2j+l n-1 

(110) = 7n-2 ^ ^!w2-2 C^'" ^/l2j;fi_i52j+l-i,n-l-p 

j=0 j=0 p=0 

n/2-1 n-2j-l 

_ \^ ^(1) ^fc,n-2i-l^fc,n-l, fe+fc 

— 7n-2 ^ ^nl2-\-'j,nl2-2^i '^i+2j ""n-l 



2j-i- 

j=l i=0 



From fl95ll we see that 



n/2-1 n-2j-l 

(111) 2^L,(A:,(n-l)/2)=7n-2 E £ ei%_,/,_^.„/,_3/,C^~^^--^C^^^^^ 

jr = l i=0 

From equation (lA6aP we have that 



(112) 



"n/2-l/2-j,n/2-3/2 ^n/2-l-i,n/2-2- 



Thus, the right-hand sides of equations f lllOp and fillip coincide. 
In the proof of the duahty (llOSbp one has to show that 

2^w(^,(^-l)/2) = Oj„,(fc,n). 

The strategy is the same as for the Laguerre ensemble and we omit the computation. □ 

We are now left with the task of computing the integrals f ll06p . When A; is a positive 
integer, we find a single sum containing k terms for each ensemble; when k is negative the 
sums go up to order of the matrix dimension. 

6.2. Incomplete Integrals — Positive Moments. We now assume that A; is a positive 
integer and focus on the Laguerre and Jacobi ensembles. 



Lemma 6.2. We have 



and 
(114) 

'2A;\ (a + 6 + 2n - 4j - 1 + 2A;)(i(a + 6 + n))(fc_j)(i(6 + n))(fc_j) 



1111111^,1,/ i — i,^; /C)1\ 



^2jJ {a + b + 2n-2j - l)(2fe+i)(l(a + n + l))(-j)(i(l + 



+ 0ln, 



where 



7L r(n/2 + fe/2-i/2) ^ ru + knb + j + k)2-^ 



r(n/2)r(6 + n - 1) ^ r(j + A; - n + l)r(6/2 + 1/2 + j) 
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and 



(116) 



2''+^T{a/2 + b/2 + j)T{a/2 + b/2 + l/2 + n/2) 
jr[ T{a + b + n + ] + k)V{] + k + l- n)V{b/2 + 1/2 + j) 

r(a/2 + 1/2 + n/2)r(& + k + j)r(j + fc) 
^ r(n/2 + 6/2)r(a/2 + l/2)r(n/2) 

Furthermore, if n > 2k, 0fe„ = 0fc„ = 0. 

Proof. We begin with the proof of formula flllSp . Equation fll06p becomes 

(117) h,{k,n) = 7.-2r/« _2 / x'^^'-'^/'C':,^^{x)e [t^'-^y^e-^/^] {x)dx. 

Jo 

The e-transform appearing in the right-hand side is the difference of the two incomplete 
Gamma functions 

/•oo 

4-0.-1 „-t. 



:ii8) 



e^^e'^dt and r(a, z) 







f-^e-'dt. 



The main idea here is to expand them in a sum of incomplete Gamma functions whose 
weight has been perturbed by a factor t''. To this end we insert 



:ii9) 



where d'j = 2^ ^ r(fe+j) • When = 1 this identity can be found in [l|, equations (6.5.21) 
and (6.5.23); the formula for general k is obtained by iteration. 
This leads to two integrals. The first one is 

k 



K,n — ''71/2-2^"- 



df+')/^x'+^^'-'e-^C'^^,{x)dx 



(120) 



(1) 



1n-2 J2 

3 = 1 

k j-\-k—\ /too 

^Y.Y. Ci^'-''''-'df^'y'x'^^^'-'e-^Ci'-\'-'r\x)d. 
,=1 i=o 



where we inserted the connection formula flM]) . Because of the orthogonality of the Laguerre 
polynomials the only contribution to the inner sum occurs at i = n — 1; furthermore, since 



max(i) = 2/c — 1, we have that 



if n > 2k. If n < 2k, inserting the appropriate 



coefficients (see appendix A) and using the duplication formula (see [3], §6) 

(121) V{2z) = t:~^'^2^'-^V{z)V {z + \) 

gives (11151) . 

The remaining non-trivial integral is 



k,n 



c 



X 2 e ^ ^n—l 



+k 

min(n— l,2fc) 



x]e 







(122a) 
(122b) 



X 2 



j=0 
min(n/2— 1,A;) 

c ^""r" 

j=0 



+fcg-./2^2fc,n-1^6+2fc_^^^^^ 



X 



x)dx 
x)dx, 
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where C = ?7^^^2-27'i-2'^fc'^^''^^/2- To obtain (I122al) we used the connection formula, while 
(]122b|) follows from the fact that the contributions to the sum with odd indices vanish due 
to the skew-orthogonality constraints (l68|l . 

Now we integrate (1122bp by parts using the identity (1791) and the formula 



6-1 



(123) -^e t'^+'^e-'/^ (x) = -x'^+'^e-^/^. 



This leads to 



dx 



min(n/2— Ijfc) n/2— j — 1 



(124) ^'^^ 

^ ' mm(n/2-l,A:) 

_ --^ \^ ^2k,n-l {i),b+2k ,b+2k 
~ ^ '^2j %,n/2-j-l"'0 

j=0 

Inserting all the relevant formulae for the orthogonality norms and connection coefficients 
gives dm]). 

We sketch the proof of flll4p as it follows a similar pattern; we only emphasise the 
differences. For the Jacobi ensemble the e-transform is expressed in terms of incomplete 
Beta functions; thus, we replace ( 11191) with the following identity: 

k 

(125) e[t\l - ty]{x) = ^d;''x'+^'(l - xy+^ + {a + b + k + l)dfe[t'+''{l - ty]{x), 

i=i 

where 

n 9fi^ - r(a + & + j + i)r(6 + i) 

^ ' ^ T{a + b + 2)Tib+l+jy 

Equation (11251) can be obtained by iteration from formulae (26.5.15) and (25.5.16) in [l[. 
Proceeding as for the Laguerre ensemble and using formula (11211) gives (I114p . □ 

6.3. Incomplete Integrals — Negative Moments. When the moments are negative, 
we focus only on the physically interesting case of the Laguerre ensemble, which leads to 
moments of the proper delay times. 

When the moments are positive the correction terms and <^^^ vanish if n > 2k. 
Now we have a similar contribution, which we shall denote (f)^^. ^ and which is not zero for 
n > 2k; however, it turns out that — )■ exponentially fast as n — )■ oo. 

Lemma 6.3. Let k be a positive integer. One has 
where 



k,n^ 



k~l 



,L _ nn/2 + b/2 - 1/2) Tjk + J + n)Tib - k - j)2^ 
^ ^ r(r2/2)r(6 + n-l)^r(6/2 + l/2-j)r(fc + j + l)' 

Furthermore, we have 

(129) = O (e-'^") , n^oo, c> 0. 
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Proof. The incomplete integral now becomes 



(130) 



X 



As for positive moments, we insert into (11301) the identity 

k-l 



(131) 



(x) 



j=0 



(6+l)/2^ 



t — 



[x)dx. 



-k-t/2 



\X 



where ci^ = 2^^^ r(fe-j) • Formula (I13ip is obtained in the same way as equation (I119p . This 
gives two integrals: the first one is 



k-l 



(132) 



= -ri%.2ln^2 ^ / S^^'^^-^-'-^C-^ CU{x)d. 



X 



k-l n~l 



_ (1) 



b-j-k-l. 



x)dx, 



j=0 i=0 

where we inserted the connection coefficients Cj for the Laguerre polynomials (16 ip . Appli- 
cation of orthogonality implies that the only contribution to the inner sum occurs at i = 0, 
yielding equation f)128p . The remaining non-trivial integral is 

1 



(-'-'^3) 2^n/2-2^n-2(lk 



X 



-k „—x 



X e 



{x)dx, 



which can be computed in the same way as the right-hand side of equation (11241) . 
Remark 6.4. The sum 

'2k + 2j - 1\ in/2 + b/2)^_k.j) 
2 J J (n/2 + l/2)(_. 



□ 



(134) 



j=0 



i) 



in equation (I127p is 0{n^^) and of subleading order compared to M^^\—k,n — 1), which 
gives the main contribution to the moments of the proper delay times for /3 = 1. However, 
it goes to zero much more slowly than the correction term 0^1^ ^. 

6.4. Gaussian Orthogonal Ensemble. The treatment of the GOE by our method is 
slightly different from the LOE and JOE. We do not ffnd a duality relation similar to 
lemma 16.11 However, the following proposition is the analogue of equation (llOOp for the 
GSE. 



k \{n/2-i- j)i^k+i/2) 



Lemma 6.5. Let k he a positive integer and suppose n is even. Then, the integral (llOSp 
is explicitly given by 

min(n/2— l,fc) min(n/2— j — j) 

(135) 0G{2k,n)= Yl E 

j=l i=0 

Proof. By changing variable of integration and using the relation fl66l) we obtain 

.00 «/2-2 



^Jv+jJ (^/2-j)(i/2) 



0G{2k,n) 



(136) 



00 



X 7n-2e «^j,n/2-2 



{Hn-i{.x)'H2j+i{.x)dx 



j=0 
oo"/2-2 



E ^ 

j=0 



fc+1/2 



7n-2e ^eV/2-2^S-i(^)^i^"(^)^^- 



,1/2 , 
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Inserting the connection formula flM]) into fll36p leads to 

n/2-2 j 

OG(2fc,n) = 7n-2 Yl 4w2-2EC^' 

(137) 

^ ^ n/2-1 



„_r, Jo 



p=0 

By using the orthogonality of the Laguerre polynomials and rearranging the indices we 
obtain 

n/2-l n/2-j-l 

(138) Oc(2t,,,)=ow.E4^2-,-,„/2-. E Cf'"''"'"'C^f"C5-.-r 

j=l j=0 

Inserting the appropriate constants from appendix A completes the proof. The coefficients 
Qk,n/2-j-i ^ those for the Laguerre polynomials; the constants 

7„_2 and e|^y2-j-i n/2-2 ^'^^ those associated to the GOE. □ 

The remaining task is the evaluation of the integral lQ{k,n) in fll06p . As previously, 
we obtain slightly different expressions depending on whether n < 2k or n > 2k. For the 
latter inequality the formula simplifies considerably. 



Lemma 6.6. Let n be an even integer. If n < 2k we have 

{2j+2i+l)j\ 



, . fc-n/2 n/2-1 /n-l\22^2!!(zl)! 
lG _ r (c., ^ _on/2-fc l^fcJ' V 2» ) {2j+2i+l)j\ 



When n > 2k we obtain 



j=0 i=0 

, (2A:)! "^^^ (n/2-z-l)!(J-/_i) 
r(n/2) ^ ^ (^■_,)!(A;-j)!4fc- 



G _.oM,^ W2 + V2-j)(,)23^-- 
(2j)!(A;-j)! 



(139) C = (2A:)!$^ 

Proof. We have 

/OO /"OO 
x2'=e-^'/2^n_i(x) / e-*'/2sgn(x - t)dtda; 
00 J —00 

where C = ''?i/2-2^"-2/2 = (2-\/7rr(n/2))^^. Now, consider the generating function 

/OO /"OO 
e^"e-"'/2Hn_i(a;) / e-*'/2sgn(x - t)rftcia;. 
-00 J —00 

Completing the square in the exponent and changing variables leads to 

(141) Mg(s) = C^e^'/M ^^.^^^)^n-i-ij / e-(''+^)'/2gg^(M - ^;)rfydM, 

j_Q ^-00 \ 3 J J -co 

where we have applied the connection formula j2^ 

n-l 



3=0 V / 
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Equation fll41p motivates us to study the function 

/OO POO 
-OO J ~oo 

In appendix B we compute fj{s) in terms of the power series 

OO OO 

(143) f,^{s) = Y,s'^+'a+^, and = 2j!v^ + J] A" 

p=j p=j+i 

where 

(UA) a+ - ^'-''V^(-^y-' and a" - ^^"^^^^(-l)^"^ 

^ ^ (2p + i)r(p-j + i)' "p'^ - (2p)r(p-j) ■ 

Thus, we can write the decomposition 

(145) Mg{s)=C{MUs)+M^{s)), 

where 

n/2-l 



(146) Mtis) = g ^^'^^(''2, ^)^"'^^-V2,(.), 

n/2-1 _ 

(147) A^c(^) = E ^^'''(2" + i)^""^'"'/2.>i(^)- 



Computing the Taylor expansions of these functions is a routine (though tedious) exercise. 
Eventually, we obtain 



(148) ^d^^ois] 



n/2-1 fn-l\ 2-J-^'(-iy 



(2A:)! '^^^"^^ (T/) 



^ 0«/2-fc y^'^)'- Sr^ Sr^ V 2i / {2j+2i+l)j\ 

nn/2) ^ ^ ik-n/2-j)l 

which vanishes if n > 2k. We also have 



■2k 



d' 

(149) ^^-^g(^) 



s=0 



{2k)l -'"^"f -^'^'^ ^ (n/2-.-l)!(J-/-i) 
r(^/2) ^ (j-2)!(A;-j)!4^-^ 



which combined with f ll35p gives the equation fl34|) . Under the assumption n > 2/c, equa- 
tion fll49p can be simplified further, leading to 

I(2kn)- ^2')' f f 4-^(./2-.-l)!(„_yj 



(2fc)!E 

j=0 



(n/2 + l/2-j)(,)23J^ 
(2j)!(^-j)! 



where in the first equality we interchanged the order of summation, while in the second 
one we have used the formula 

(150) t ' 



ik - j)\u - py. {k-p)\- 



□ 
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Appendix A. The Coefficients /i„, and e 



(4) 



The orthogonality normahsations hj defined in fl551) for the classical orthogonal polyno- 
mials are tabulated in any standard reference on special functions jH |52|. However, in this 
article we work with monic polynomials. This is not the standard normalization found in 
the literature. Therefore, for the reader's convenience, we report the coefficients hn, e 



Cn^ that we use throughout this paper. 
We have 



(1) 

n , 



(Al) 



ua,b ^ r(a+i+l)r(6+i+l)r(j+l)r(a+b+j+l) 
r{a+b+2j+l)r{a+b+2j+2) ' 

/l5 



r(j + i)r(6 + j + i). 



Jacobi, 

Laguerre, 

Hermite. 



Given these normalisations, recall that 



(A2) 
where 



(A3) 



hnjn 



1, Hermite, 
|, Laguerre, 
|(2?2 + a + 6 + 2), Jacobi. 



In §1] for the orthogonal and symplectic ensembles we introduced the quantities 



(A4) 



h 



2n+2 



n 



C2i 



%•+! C2i+l 



.(4) 



n-1 

't'2n+l TT C2j+1 



C2r. 



n 



C2i 



n 

j=0 



C2jh2j+2 
C2j+lh2j 



We have the following explicit formulae: 
JOE and JSE, 



(A5a) 



(A5b) 



(A5c) 



,(4),a,f> 

"j,n 



X 



16"~-'2r(n + 3/2)r(n + 3/2 + a/2) 
r(j + 3/2)r(j + 3/2 + a/2)r(j + 3/2 + 6/2) 
Tin + 3/2 + a/2 + 6/2)r(4j + 4 + a + 6) 
r(j + 3/2 + a/2 + 6/2)r(4n + 5 + a + 6) 
16"-J'2r(n + l)V{n + a/2 + 1) 
r(j + l)r(j + l + a/2) 
T{n + 6/2 + l)T{n + a/2 + 6/2 + l)r(4j + a + 6 + 2) 



X 



r(j + 1 + 6/2)r(j + a/2 + 6/2 + l)V{An + a + 6 + 3) 
r(a/2 + 6/2 + 3/2 + n)r(6/2 + 3/2 + n)r(a/2 + 3/2 + n) 
2-4n-4-a-b7^r(^Q + 6 + 4n + 5) 

r(a/2 + 6/2 + l)r(n + 3/2) 



X 



r(6/2 + l/2)r(a/2 + l/2) 
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LOE and LSE, 




(A6a) 




(A6b) 




(A6c) 




GOE and GSE, 




(A7) efl - 


V{n + . 



4"-^2r(n + 3/2)r(n + 3/2 + 6/2) 

r(j + 3/2)r(j + 3/2 + 6/2) 
4"-^'2r(n + l)r(n + h/2 + 1) 

r(j + i)r(j + 1 + 6/2) 

^^, r(n + 3/2)r(n + 6/2 + 3/2) 
V^r(6/2 + 1/2) 



vr 



Appendix B. The Generating Function fj{s) 
In the proof of lemma 16.61 we needed to study the generating function 



(Bl) 



^(''+s)V2ggj^(^y - v)dvdu. 



The signed integral in the equation fIBip is closely related to the error function, for which 
the k^^ derivative can be expressed in terms of the monic Hermite polynomial T-Lk-iiu)- 
Differentiating under the integral, we find 



(B2) 



1 d 



2k 



{2k)\ds^'' 



f2jis) 



s=0 



{2k)\ 



V2 e''''n2j{u)n2k-iiu/V2)du = 0, 



which follows from the oddness of the integrand. 
For the odd derivatives we get 



■2k+l 



{2k + 1)\ ds^''+ 
(B3a) 

(B3b) 



-Ms) 



2fc+i 



poo 

— / e~^"n2jiu)n2kiu/V2)du 

e-'^u-^'^ q^'^{u)C'^^'\u/2)du 

-^u-^'^C-:''\u)CT''\u)du 



(2/c + l)!.o 
^(k-l/2\k\ 2 

1=0 



k-i J i\ (2A; + 1)! Jq 



^^^"^ {2k + l){k-j)\- 

The expression in line (lB3aP was obtained by means of the relations f l66p . while (lB3bp 
follows from applying the connection formula 



(B4) 



£|(x/2) = 2-'$:(tt*)|(-i)"':?W- 



i=0 



The last line (IB3cp then follows from the orthogonality ( l38p . Thus, we have 
(B5) 



(2p+i)r(p-j + i)- 
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In a similar fashion we find 

(2p)r(p-j) ' 



(B6) f,,,M = 2,!0F + Y: Jv ~- \^ 



p=i+i 

which required the double integral 

/oo /"OO 
-00 J —00 
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